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Positive parity meson states -Dso(2317)''' and -Dsi(2460)''" have masses slightly below the DK 
threshold. Both states can strongly decay only into isospin violating decays Dsi(2460)’^ —>■ Dfmr, 
P>si(2460)''' —>■ and Z)so(2317)’^ —>■ The tt states have rather small energies, which 

makes these decays appropriate to study within heavy meson chiral perturbation theory and calcu¬ 
late loop contributions. The Dsi(2460)''" —>■ decays occur only at the loop level. Amplitude 

is a result of chiral loop contributions, which then have to be finite. However, in the case of 
i7si(2460)''' —>■ and D*q{2317)'^ —>■ decays, there is a tree-level contribution. We find 

that chiral loop contributions might be important in both cases. The calculated amplitudes are sen¬ 
sitive on the coupling constant describing the interaction of positive and negative parity heavy meson 
multiplets with the light pseudoscalars. The counterterms contributions are also present in the am¬ 
plitudes Z)si(2460)^ —>■ D*+7r° and Z)*o(2317)^ —>■ We explore an experimentally known ra¬ 

tio of the decay widths for these two decay modes to estimate the size of counterterms contributions. 
We determine decay widths for both decay modes to be r(Dsi(2460)’^ —>■ ~ 0.25 keV and 

r(D,i(2460)+ D+7r°7r°) ~ 0.15 keV. 


* Electronic address: svjetlana.fajfer@ijs.si 
^ anita.prapotnik@um.si 



2 


I. INTRODUCTION 


It has been more than a decade since Z?*q(2317)“'" and ZIsi(2460)''" mesons were observed by the BABAR and 
CLEO collaborations, respectively. Their existence was confirmed by several experiments [l]-[5]. The measurement 
of D*q{2317)'^ and I?si(2460)''' quantum numbers, as well as decay widths and decay rates, has continued since 
that time [l]-[9]. The experimental data support the interpretation of the ZI*g(2317)^ meson as a positive parity 
scalar (J^ = O’*"), while ZIsi(2460)''" appears to be positive parity axial vector (J^ = I"*"). Both states behave as cs 
systems, although their masses turned out to be 100 MeV smaller than expected by calculation based on the quark 
models (for a review see [10]). Many proposals have suggested that these states are tetraquarks or DK molecules 
[10]-[14]. There is also a suggestion that these states can be a result of the mixing of sc and scqq states (see e.g. [15]). 
Recently, lattice studies indicated that II*q(2317)+ is a combined state of sc and DK molecules [16-19]. However, 
multiple lattice volumes will be needed [18] to resolve the structure of Zl*g(2317)+. On the other hand, the total 
decay widths of II*q(2317)+ and IIsi(2460)+ are still unmeasured. The upper bound on the total decay width of 
the Zl*g(2317)''' meson is 3.5MeV at 95% confidence level, while the upper bound on the total decay width of the 
Ilsi(2460)+ meson is 3.8MeV at 95% confidence level [7]. Nevertheless, some branching rations for IIsi(2460)''' 
decays were determined. The masses of II|q(2317)+ and ZIsi(2460)+ mesons are slightly below the threshold for the 
decay into a D and a K meson. Therefore, only strong isospin violating decays of ZI*g(2317)+ and ZIsi(2460)+ are 
kinematically allowed, as well as radiative decays of II|q(2317)+ and IIsi(2460)+. The decay II*g(2317) —?► D'^tt^ was 
observed by the BABAR collaboration [1], but branching ratio or partial decay width was not determined. However, 
the branching ratio of I?si(2460)''' —)■ was found to be (48 ± 11)% [7, 20]. Thr BABAR Collaboration also 

observed IIsi(2460)+ —>■ D+tt+tt” [6], with the branching ratio of (4.3 ± 1.3) % [20]. There is no measurement of the 
branching fraction for the decay mode Zlsi(2460)''' —>■ yet. In the case that T)*g(2317)’^ and ZIsi(2460)'*' are 

only cs mesons, with quantum numbers 0+ and 1+, the decay of II*g(2317)+ into is possible, while the decays 

of II*g(2317)+ into II*+7r° and H+tttt cannot occur. Similarly, the decay of Zlsi(2460)''" into D'^tt^ is not possible, 
while the decays of Dsi{2AQQ)^ into ZI*+7r° and are allowed. Therefore, it seems that existing experimental 

results favor the treatment of both states II*g(2317)’^ and I?si(2460)“'" as a cs states with the positive parity. 

The T>si(2460)+ —)■ ZI*+7r° and T>*g(2317)+ —>■ Df 7r° decays were considered within variety of approaches, [21]-[33]. 
The quark models were mostly exploited for the decays II*g(2317)+ —)■ Dfr] —)• and Zlsi(2460)''" —>■ Dl~^r] —)■ 

D*'^Tr^, assuming isospin symmetry violation via rj — tt mixing. In Ref. [22] states ZIsi(2460)“*' and Zl*g(2317)^ are 
treated in the DK- molecule picture. Most of these studies predict the partial decay widths of Zlsi(2460)''" — 
and II|g(2317)+ —7r° decay modes to be in the range (10 — 30) keV, although higher values (about 100 keV) were 
suggested in Ref. [24] too. 

The three-body decays T>si(2460)+ —)■ D+tt+tt” and ZIsi(2460)+ —)■ ZI+7r°7r° were studied in Refs. [34, 35]. The 
authors of [34] assumed that the decay occurs trough the intermediate a fields, which then convert to pions through the 
ss component of a fields mixing with the uu and dd components. Two relevant cr states (erg and erg) were considered 
with masses set to 1 GeV and 1.5 GeV. In Ref. [34] high sensitivity of the amplitude on the mass of the lighter sigma 
meson state was found out and the variation of its mass in the range 0.8 GeV and 1.2 GeV can change the predicted 
result for an order of magnitude. In [35], the a states were replaced by the scalar /g(980) state. However, the decay 
width was estimated to be r(ZIsi(2460)+ —>■ D+tt+tt”) < 25keV. 

In this paper, we determine chiral loop contributions to the isospin violating decay amplitudes of ZI*g(2317)'*' and 
Ilsi(2460)^ mesons. For two-body decays, there is a tree-level contribution to decay amplitude arising from the rj-ir 
mixing. However, chiral loops even in these cases, might give significant contributions. This was indicated already in 
Ref. [30], where some of the loop contributions to the decay amplitude were determined. In our analysis, we rely on 
the heavy meson chiral perturbation theory (HMyPT). Within HMyPT the Zl*g(2317)+ and IIsi(2460)“'' states have 
quantum numbers of cs. The use of HMyPT in both decay modes is fully justified by the fact that the pions in the 
final set have rather small energies. The loop contributions within this framework arise from the light pseudoscalar 
meson exchanges. The light resonances, as light vector mesons (p, K*) in the amplitudes at tree level, give the 
contributions of the same order in the chiral counting [36] as light pseudoscalar meson loops [37]. In comparison with 
the approach of [30], we find that there are additional Feynman diagrams leading to the relevant contribution to the 
two body decay amplitudes. The three body decay amplitude within this framework arises from chiral loops. The 
energy release in both two-body and three-body decays is very small. Both negative and positive parity intermediate 
D states are taken into account within this framework [38-40]. 

The heavy meson Lagrangian formalism will be introduced in Sec. H. In Sec. HI, the analysis of the Z?si(2460)“'' —>■ 
D+tttt decay channels will be presented, while in Sec. IV we will calculate the decay width of the two-body 
Zlsi(2460)+ —>■ ZI*+7r° and ZI*g(2317)+ —>■ Dfn^ decay modes while a short conclusion will be given in Sec. V. 
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II. FRAMEWORK 

The framework of heavy meson chiral perturbation theory combines the heavy quark effective theory with the chiral 
perturbation theory [41, 42]. Heavy quark effective theory is used to describe mesons composed of one heavy quark 
(c or b) and one light quark (u, d or s) [41, 42]. In such mesons, the heavy quark moves almost on shell with the 
velocity and the momenta of the heavy meson can be written as = mv^ + where m is a heavy meson mass 
and is of the order of Aqcd and much smaller then mv^. In a limit, when the mass of the heavy quark becomes 
infinite, pseudoscalar and vector meson states become degenerate, as well as scalar and axial vector meson states. 
The negative parity states are described by the held H, while the positive parity states are described by the held S: 

H=^{l + v-t)[P;r-P75], S=^{l + vj)[D;rj,-D], ( 1 ) 

where P* and P annihilate the vector and pseudoscalar mesons, respectively, while and D annihilate the axial- 
vector and scalar mesons, respectively. 

Within chiral perturbation theory, the light pseudoscalar mesons are accommodated into the octet S = P 

with 


/tt^/v^-I- 778/-\/6 7r+ Ar+ \ 

n = I 7r“ -7r°/v^-|-rys/Vb K° (2) 

V K- K° -2r]s/V6j 

and / ~ 120 MeV at one loop [43]. The leading order of the HM^PT Lagrangian, that describes the interaction of 
heavy and light mesons, can be written as 

C = -Tr[Ha{iv ■ Vab - Sab^H)Hb\ + gTr[HbHa'r ■ Aablb] 


+ Tr[Sa{iv ■ Vab) - 5ab^s)Sb] + gTrlShSa-^ ■ Aablb] + hTr[HbSa'^ ■ Aablb] , (3) 

where = Sabd^ — Vab ^ heavy meson covariant derivative, is the light meson vector 

current, and Afj, = is the light meson axial current. A trace is taken over spin matrices and 

repeated light quark flavor indices. All terms in (3) are of the order 0{p) in the chiral power counting (see, e.g.,[39]). 
As in [39] we assign for Ash — As — Ah ^ 0{p) in order to maintain a well-behaved chiral expansion. Light mesons 
are described by the Lagrangian [41, 42], which is of the order 0{p^) in the chiral expansion: 

>Co = + !^TT[mp: + Em«], (4) 


with Ao = ml! (m„ -I- mA) = 
man rules, as given in [39]. 
written in the forms 


(m^+ — m\ )/(m„ -I- md) = — m^l2)lms . The above Lagrangians lead to Feyn- 

The scalar [pseudoscalar) and vector (axial-vector) heavy meson propagators can be 


i 

2{k ■ V — Ai) 


and 


2{k ■ V — Ai) 


(5) 


respectively, where in the propagator represents the residual mass of the corresponding field. Residual masses 
are responsible for mass splitting of heavy meson states. The difference As — Ah splits the masses of positive and 
negative parity states. In addition, we also have to take into account mass splitting between Ds and D states as 
well as mass splitting between vector (axial-vector) and pseudoscalar (scalar) fields. These splittings arise due to the 
heavy meson Lagrangian correction of the order 0{mq) (m^ stands for the mass of light quarks). To account for all 
of the above mass splittings, we will follow the approach of [38, 39] and set the values of to the experimentally 
measured mass differences between D meson states. We use the mass of the initial particle as a reference value, so all 
mass differences are defined as mass differences between the relevant D meson and initial state [38]. 

The coupling constants g, h, and g were already discussed by several authors and determined by several methods: 
the QCD sum rules [44]-[48], the lattice QCD [49]-[55], and the extraction from the experimental data [38, 39, 56, 57]. 
We will use recent results from the lattice QCD: g = 0.54(3)(l4) [49], g = —0.122(8)(6), and h = —0.84(3)(2) [55]. 
The values of h and g were determined for the B meson sector, so l/rric corrections can make a slight difference in 
stated values. Although HMy^PT relies on the expansion in the light pseudoscalar momentum and l/rric expansion, 
we do not consider l/rric corrections for at least two reasons: first, the number of additional terms in the Lagrangian 
becomes huge and impossible to estimate and second, lattice studies indicate that these contributions are rather small 
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[53]. The authors of [55] made an estimation of h for the D meson sector and found that h ~ 0.74(8), which is still 
inside the error bars ol h = —0.84(3)(2). 

In order to absorb divergences coming from loop integrals, one needs to include counterterms in the Lagrangian. 
Following [38, 39], they can be written as 

Cct = \l[HbHa{m^q)ba] + Aj [.ffa^a (w^ )bh] - A~i [;5&5a(TO|)ha] “ A'^^ [;5a5a (to^ )bh] + 

^^^Tr[(i75'7^y4'"75)ab(TO«)f,a] + ^^^^^Tr[{HS-i^A^^-ib)aa{ml)bb] + 

^^rr[iJa5'a7MAlbc75(TO^)cb] + J^^^Tr[H^Salf,A'i^-ib{ml)ab] + 


j^^Tr[HaSbiv ■ Vbcl^,A^alz] + J^^Tr[HaSbi-^^, ■ + h.c. + ... , ( 6 ) 

where = {^niq^ — ^^mq^^)/2 and Dbc^ca = ^'^^ba + b“^^]ba- At the given scale, the finite part of Kg can be 
absorbed into the definition of h. Parameters Ag and A)^ can be absorbed into the definition of heavy meson masses 
by phase redefinition of H and S, while Ai and Ai split the masses of SU(3) flavor triplets of Ha and Sa [38, 39]. 
Therefore, only contributions proportional to k[, Kg, Kg, S 2 , and S'^ will be explicitly included in the amplitudes. 


III. THE L>^i( 2460 )+ ^ D+tt+tt' AND Dsi( 2460 )+ ^ D+ 7 r° 7 r° DECAY MODES 


First, we will consider the chiral loop contributions to the Z3si(2460)“*' —>■ D+tt+tt decay rate. The decay width 
for the Dsi(2460)+ —)■ D+tt+tt”, averaged over the Dsi(2460)’'' polarizations, can be written as: 


dr = 


1 1 
(27r)3 32Mf 


|A4|^dmi2dTO23 ^ 


(7) 


where Mi is the mass of initial particle Dsi(2460)3', dmfg = {p+ +P-)^, and dm^g = (p_ + q)^. Here p_, and p+ 
denote the momenta of 'k'^ and 7 r“ respectively, while q is the momentum of D'^. The decay amplitude Ai, in general, 
can be written in the form 


Ai = e ■ {Aip+ + A 2 P- + A^q). 

Within heavy quark limit = MiV^, q^ = MfV^, and p^ + p^ = Amv^, where Am = (Mi — Mf) and Mf is the 

mass of Df meson. As e • u = 0, the amplitude simplifies to 

Ai = Ae ■ (p+ — P-) = Ae ■ Ap , 

where A can be calculated from the diagrams presented in Fig. 1 . Only diagrams giving a nonzero contribution 
are shown. Note, that there are no diagrams with the 77 meson in the loop on Fig. 1, as they all give a vanishing 
contribution. 

The diagrams in Fig. 2 also give a nonzero contribution to the amplitude. Since they are next-to-leading order 
in HMyPT their contributions can be neglected. Note, also, that by taking = Miv^ and = MfV^, all 
scalar products of the momenta become independent of the phase space parameters: ■ v = p- ■ v = AM/2 and 

P- ■ p+ = AM/2 — m?. Therefore, the spin averaging of the amplitude Ai is constant on the whole phase space 
region, implying that the calculation of the amplitude and the integration over the phase space in (7) can be done 
independently. 

Now, we proceed to the calculation of the amplitude A. Using Feynman rules derived from (3), we obtain 


where 


A = 


hJM,Mf , 

207 ^ 2^4 -f 02 -|- 61 -|- 62 -|- Cl -f C2) , 


Oi = I {Bi{Ap,,mKo) - Bi{Ap,,fnK+)) , 


( 8 ) 


( 9 ) 
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Figure 1. Non-zero contributions to Dsi (2460)"'" —>■ D^-k^-k decay amplitude. 
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Figure 2. Feynman diagrams contributing to Dai(2460)"'" —>■ D+tt+tt decay amplitude at higher order and which 
included in our calculations. 


02 = ^ (.Bi(A_d -f Am, rriKo) — Bi{Ad + Am,TO/ f+)) , 


^1 = ((i?2(Ap*, Ap -b AmI‘^, rriKo) — AmI‘^ • (Ap., Ap -b AmI‘ 2, rriKo)) 


~ (-B2(Ap., Ap -b Am/ 2, mp:+) — AmI 2 ■ B[{Ap,, Ap + Am/ 2, TOp-+))) , 


^2 = 2g ((i32(Ap. -b Am 12, Ap + Am, wp-o) -b Am 12 ■ (Ap. -b Am/ 2, Ap -b Am, mKo)) 


are not 


( 10 ) 


( 11 ) 


(i? 2 (Ap. -b Am/ 2, Ap + Am, rnji+) -b Am 12 ■ B[{Ad* , Ap -b Am, tnp:+))) 


( 12 ) 
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Figure 3. Tree-level contribution to -Dsi(2460)^ —>■ -k^ and Dso(2317)''' —>■ D decay modes. 


Cl = —2g (^(^Boo{mKo) + Ap.C'oo(Ap- , m^o)) — (i?oo(™if+) + Ap.C'oo(Ap.,mp-+))) 


C 2 — —2g ((i3oo(wpio) -f ApC'oo(Ap, TOp-o)) — (i?oo (top:+)-I-A pC'oo(Ap,TOp:+))) . (13) 

Here, Bi, B 2 , Bqq and Coo are the Passarino-Veltman loop integrals defined in Appendix A. Note that these loop 
integrals, as well as the parts of the amplitude oi, hi, and Ci, are divergent in the H > 4 limit. Nevertheless, by 
summing up all contributions, divergences cancel out, so both A and M are finite. As I?si(2460)+ —>■ H+tt+tt” 
does not have any tree-level contributions from the heavy meson Lagrangian, this was expected. Another interesting 
feature is, that the amplitude A4 is nonzero due to the mass difference of the K'^ and mesons. Namely, if we put 
TOp-+ = mpTO, the amplitude would vanish. Finally, the decay width coming from these amplitudes is 

r(i:)si(2460)+ ^ H+TT+TT-) = 0.25(4)(7)(t4)keV. 

The first error comes from the uncertainty in the coupling constant h, the second from uncertainty in the coupling 
constant g, and the last from the uncertainty in the mass of Zlsi(2460)“'" meson. Uncertainties in the coupling constant 
g and other meson masses are relatively small and therefore can be safely neglected. This result implies that the total 
decay width of Z?si(2460)''" is found to be between 2keV and 13keV. 

The decay amplitude for Dsi{24:60)~^ —>■ Z?+7r°7r° can be easily found by replacing the tt+tt” state by in the 
final state. Note that one has to include the factor 1/2 in the decay mode due to two identical bosons in the final 
state. This yields 


r(i:)^i(2460)+ ^ T>+7r°7r°) = 0.15(4) keV. 


IV. THE L>„i(2460)+ ^ Dt+iv^ AND D/o(2317)+ ^ D+7r° DECAY MODES 


The tree-level contribution to Ilsi(2460)+ —>■ Il*+7r° and Z1 /q( 2317)+ —decay amplitudes results from 
the T] — TT mixing. In this scenario, the decays proceed through the channels Z?si(2460)“'" —?► 0*^7] and 

Zl*g(2317)+ —>■ D'^t] —)■ as presented in Fig. 3. The ry —7r° mixing can be approached by the mixing Lagrangian 

[58, 59], 


r = 

^rj — no 


ml(mu - mg) 
V^(m„ -I- rud) 




(14) 


which comes from the second term in (4). The decay width for Dsi{2460)~^ —)■ coming from this mixing 

tree-level amplitude is 


r(Z4,i(2460)+ ^ D*+n°) = = 16keV, (15) 

where and Et^ are the momenta and energy of the outgoing pion, while mixing parameter 5mix is defined as in 
[58, 59] 


. _ _ mu - md 

2V2ms - {mu + md)/2 ' 

Note that decay widths for T>si(2460)+ —)■ and D*q(2317)+ —>■ 

difference in Df masses. 


(16) 


decay modes differ only for a small 









7 



K°(K') 



K^orK' K°arK^ 

Figure 4. Chiral corrections to the -Dsi(2460)^ —>■ decay mode. 
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Figure 5. Chiral corrections to the DIq{2317)^ —>■ decay mode. 
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Figure 6. Chiral corrections to the D meson wave functions. 
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h 


h 


Figure 7. Dependence of _Dsi(2460)^ —> _DJ"'"7r° and -D3o(2317)''" —>■ -Dj"7r° decay widths on input parameters at Act = 0. 
Dashed lines present the uncertainty on the decay width due to the uncertainty in g and the mass of the initial meson. 


where Zf and Zi denote wave function renormalization of the initial and final heavy meson states, 




^ ^ _ l dTl,j{vp) 


2 dv ■ p 


on mass shell 


while Zy represents the vertex corrections 


Zy = l- 


t{v Pi,v Pf^k"^) 


ro{v ■ pi,v pf,k‘^) 


on mass s. 


hell 


(18) 


(19) 


Here F is the vertex amplitude calculated from the Feynman diagrams presented in Figs. 4 and 5, while Fg is the 
vertex amplitude resulting from the tree-level Feynman diagrams (see Fig. 3). Similarly, n(z; ■ p) is the heavy meson 
self-energy arising from the Feynman diagrams in Fig. 6. 

The vertex corrections come from the Feynman diagrams presented in Figs. 4 and 5 and can be summarized into 
the expression 

Zy = l- - ^{A^{mK+) + AiirriKo)) + -^=^—(A: (TOic+) - A(miyo)) -h At) , (20) 

where = 0.11 includes corrections to the ry — tt mixing angle beyond tree level [38, 43] and the functions Ai are 
given in Appendix B. 

The isospin violating nature of both decays are manifested either by proportionality of the amplitude to the mixing 
parameter Smixj or by vanishing of the amplitude in the case of the isospin limit m^o = mx+- 
Finite parts of counterterms, are included in amplitude as Af 


Act — 


327r2/^ 




4) + 


2 2 I '42(r 


i'K+ 


i-KO) 


4 + “^(4 + 4) 


( 21 ) 


The values of the finite parts of counterterms, of course, depend on the renormalization scheme. We use dimensional 
regularization in the renormalization scheme in which the divergence 2/e contains the constant — 7£;-|-ln47r-|-l, coming 
from the loop integrals. This has to be taken into account when discussing the numerical value of the Act term. The 
wave function renormalization terms Zf and A arise from sunrise diagrams presented in Fig. 6 

Zj = 1- Bj{mK+) - BjirriKo) - ^B'j{m^). (22) 

The functions Bj are listed in Appendix B. 

The calculated numerical values of the decay widths of Zlsi(2460)+ —>■ Zl*+7r° and i4*g(2317)+ —>■ are very 

sensitive on the value of h and can be modified by two orders of magnitude when h varies form —0.65 to —0.9. The 
wave function corrections Zij are source of that sensitivity. The decay widths are also moderately sensitive on the 
mass values of final and initial states and vary on coupling constant g. The dependence on the decay width on g and 
the masses of the final and initial mesons are presented in Fig. 7, while the decay widths of Ilsi(2460)’'' —)■ 
and Zl/o(2317)+ —>■ are presented in Fig. 8, for a range of values for counterterms At ■ 
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h 


Figure 8. Left: Dependence of Dsi(2460)^ —>■ and Dao(2317)''" —^ decay widths on the decay constant h for hve 

different values of counterterm amplitude. From thickest to thinnest line, Act ~ 0,—1,1,—3, 3. The Dsi{2AQQ)'^ — >■ L>*+7r° 
decay width is presented by the solid line, while the L)so(2317)^ —>■ decay width is presented by the dashed line. Right: 

Dependence of Dsi(2460)^ —^ and D*o(2317)^ —> decay widths on the value of counterterm amplitude for three 

different values of coupling constant h. From thickest to thinnest line, h = —0.89,—0.84,—0.79. The solid line denotes the 
decay width of Dsi(2460)''' —^ while the dashed line is used for the decay width of Dso(2317)''' —>■ Dj"7r°. 



Figure 9. Determination of the allowed regions for the counterterm size in the amplitude Act- Horizontal lines present the 


experimental values of 
calculation for the ratio 


Sfl(Dal(2460) + -»D* + 7r°) 


^ with the one sigma error band. The ”U” line presents the result of our 


BR{D,i{2A60) + 

of h in the range -0.79 and 0.84. 


SH(Dai(2460) + -!'r3a + 7r+7r-) 

Bfl{Dsi{2460) with the band region presenting the uncertainty coming from the variation 


-) 


Since the numerical results are very sensitive on the value of coupling constant h and due to the unknown final parts 
of counterterms, we cannot make a definite prediction for the partial decay widths for both decay modes. Nevertheless, 
by using experimentally measured ratio of these rates. 


Bi?(£)si(2460)+ ^ Dl+Tr°) 
BR{Dsi{2m)+ Dt+TT+T:-) 


0.09 ±0.02, 


(23) 


we can shed more light on the value of counterterm amplitude Act- Our result given in Fig. 9, indicates that Act can 
be accommodated within the ranges (-4.8;-3.2) and (-1.3;-0.1). Here, we only considered one tr experimental error 
on the decay widths ratio and the uncertainty in the theoretical prediction, by varying h between —0.79 and —0.89. 
Within these bounds, we obtain for the Dsi(2460)+ —)■ D*+7r° and D*q( 2317)+ —)■ decay widths: 

r{D,i{2460)+ D*+Tr°) = 2.7- 3.4keV, (24) 

r(i:)*o(2317)+ ^ i:»+7r°) = 2.4 - 4 . 7 keV . (25) 


One can notice a slight difference between r(Dsi(2460)'*' —>■ D*~^Tr°) and r(I?*Q(2317)’'‘ —not present at 
the first order which arises from the r] — mixing. This difference is a result of loop corrections, which were not 
considered in previous calculations [44]. 
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V. DISCUSSION AND CONCLUSIONS 

Within a HM^PT framework, we determine loop contributions to the strong isospin violating decay amplitudes 
for Z?si(2460)+ and Zl*g(2317)''". We have assumed that both states carry quantum numbers of the cs states. Since 
three-body decays of these states are forbidden at tree level, we calculate contributions to the decay amplitude at 
the loop level, which is consequently finite. Contrary to three-body decay amplitudes described by the finite loop 
contributions, two-body isospin violating decay amplitudes receive contributions at tree level, induced by the rj — 
mixing. Therefore, the chiral loop contribution is not hnite and in order to regularize it, one has to introduce 
counterterms. Our estimate of the size of counterterm in the amplitude Act relies on the result we derive for the 
Ilsi(2460)+ —>■ U+TT+TT” decay width. 

Note that within chiral perturbation theory, only the light pseudoscalar mesons are present in the loops. Contribu¬ 
tions of the light resonances with the spin = 1“, 0+, I'*' are already accounted by the chiral loop contributions of 
light pseudoscalar mesons [60, 61]. Additional inclusion of light resonances with = 1“, 0+, I’*' in HMyPT is not 
consistent with the original framework. Nevertheless, one can roughly estimate contribution of light vector mesons 
in the loop. For example, the contribution of K* vector mesons in loops as seen in [12] would correspond to two 
loop effects within our framework. Their contribution to the decay rates is therefore suppressed within HMyPT, in 
contrast with the K* loop contribution in the D*K molecule picture presented in [12]. 

If one assumes as in Ref. [34] that heavy mesons interact with spin = 0+ light resonances as cr (or /o(500) 
as in [20]) at tree level, then Z?si(2460)+ —D+tt+tt” decay can proceed through Zlsi(2460)''" — Df a D+tt+tt”. 
However, their contribution is proportional to the e • v which is then equal to 0 in HMyPT. On the other hand, one 
might think that since the mass of cr is close to 500 MeV, that state can be important in Z)si(2460)+ —>■ D+tt+tt”. 
This would be the case if a contained the ss component. Recent lattice QCD study [62] does not support such idea. 
In principle, higher order terms in the chiral expansion might contain terms which describe interactions of heavy 
meson states with = O'*' light resonances. At the same time that will mean that some of these fo states should 
contain ss contribution. Unfortunately, the structure of light positive parity scalar mesons is not known yet and 
reliable consideration of this contribution is not possible at present. 

Better understanding of the structure of Zlsi(2460)+ and I1 *q( 2317)+, as well as light scalar mesons, might shed 
more light on the decay mechanism for two- and three-body strong isospin violating decays of i7si(2460)“'' and 
^:o(2317)+. 
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Appendix A: Loop integrals 

By employing dimensional regularization, in the renormalization scheme with 6 = -I- ln47r -|- 1 = 0, we 

have 


Ao{m) = 




nr‘‘ J — m‘‘ + le) \ fj,‘‘ ' 


Bo{p, 171,171) = -- 


d^k 


iiT^ J {k"^ — m? + ie){{k + p)"^ — rri^ + ie) 
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1 


Booip, 717, m) = [Ao(m) -H (2m -p /2)Bo(p, w, m)] 


which in Z? —>■ 4 limit gives 


Booip,m,m) = - [Ao{m) + (2m^ - p^/2)Bo{p,m,m) + 2m^ -p^/3], 
o 
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Booim) = BooiAMV,m,m). 
Loop integrals with one light meson propagator are 

f d^k 


Bo (A, to) = 


ZTT^ J {k‘^ — m‘^ie){v • k — Aie) 


-2A 




+ 0(B-4), 


with 


Fil/x) = 


j — 1 ln(a; + — 1 + ie); |a:| > 1, 


Vl-x2(f-tan 1 , 


Ixl < 1. 




d^k 


(k^ — to2 + ie){v • k — A + ie) 


— = Bi{A,m)vf" , 


Bi(A,to) = 


(27r/i) 


4-D 


k■Vd^k 


J (k"^ — + ie){v • k — A + ie) 


— = Ao{m) + ABo(A, to) . 


B^''{A,m) = 


(27r/r) 


4-D 


d^k 


J (k^^ — m‘^ie){v • k — Aie) 


— = Boo{A,m)g>"'" + Bii{A,m)v^v'' 


■Boo(A,to) = ^ ^ [(to^ - A2)Bo(A,to) - AaIo(to)] , 


which in the D —>■ 4 gives 


_Boo(A,m) = -[(m^ — A^)BQ{A,m) — AAo(m) + 2A/3(3r?i"^ — 2A^)], 

o 


Bii(A,to) = [{DA^ - to^)Bo(A,to) + L>AAo(to)] , 


which in B —> 4 gives 


Bii(A, to) = -[(4A2 — to^)Bo(A,to) + 4AAo(to) — 2A/3(3to^ — 2A2)] , 

O 


B2(A,to) = Boo(A,to) + Bii(A,to) , 




d^k 


(fc2 — 'm?){v ■ k — Ai)(r) • k — A 2 ) Ai — A 2 


[Bo(Ai,to) - Bo(A2,to)] , 


B'''(Ai,A2,to) = 


(27r^) 


4-D 


k>^ d^k 


iir^ J {k"^ — Tv?){v ■ k — Ai){v ■ k — A 2 ) 


= B[{Ai, A2,m)v^ , 


Bi(Ai, A2 ,to) = 


{2tth) 


4-D 


k • Vd^k 


ZTT^ J {k‘^ — m‘^){v • k — Ai){v • k — A 2 ) 


= _Bo(A2,m) + AiBq(Ai, A2,m), 
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B2(Ai, A2,to) 


(27r/x)4-^ 


(k ■ vY d^k 

(fc2 — 'm?){v ■ k — Ai)(t) • k — A2) 


Ao{m) + (Ai + A2)Bo(A2,m) + AiBq(Ai, A2,m), 


B^"'(Ai,A2,m) = 


(27r/i) 


4-D 


d^k 


iTT^ J (k"^ — m'^){v ■ k — Ai){v ■ k — A2) 


S' 0 (Ai, A2, m)g^‘' + S; 1 (Ai, A2, , 

BooiAi,A2,m) = ^ [m^So(Ai, A2 ,to) - AiSj(Ai, A 2 , m) - Si(A2,m)] , 
which in S —)■ 4 gives 

i[m^So(Ai, A2,m) - AiSj(Ai, A 2 ,m) - Si(A2,m) + 2/3(3m^ - 2(Ai + A 2 + A 1 A 2 ))], 


Sii(Ai, A2 ,to) = jj [-to^So(Ai, A 2 , m) + SAiSj(Ai, A 2 , m) + SSi(A2,to)] , 

which in S —> 4 gives 

to^Sq(Ai, A2, m) + 4 AiS^(Ai, A2, m) + 4Si(A2, m) — 2 / 3 ( 3 to^ — 2 (A^ + A2 + A1A2))], 


Loop integrals with two light meson propagators are 

(27r/r)4-^ f 


A, mi, m 2 ) = 


d^k 


ITT^ J (fc2 _ 171,2 _|_ — p)2 _ 1772 _|_ ('ll . 


Cl (p, A,mi,m2)n'^, 


Clip, A, mi, m 2 ) 


ITT^ 


k■Vd^k 

(k'^ — m\+ ie)i{k — p)^ — m^ + ie)(n • fc — A + le) 


So(p,mi,m2) + AC'o(p, A,mi,m2), 


C'^''(p, A,mi,m2) 


(27rp)4-^ 


d^k 

(fc^ — mf + *e)((fc — p)2 — m| + ie)(r' • fc — A + ie) 


C'oob, A,mi,m2)p^'' + C'ii(p, A,mi, m2)?;'‘w'', 


C'oo(A,m) = Cooi-AMV,A,m,m) 


D -I [-^o(-Am + A,m) - (Am /2 + A)So(A„n, m, m)+ 
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(to^ - A^)Co{Amv, a, m, m)], 


which in the —>■ 4 gives 


Coo(A,to) = Coo(-AM'y, A,m,m) = -[So(-Am + A,m) - (Am/2 + A)Bo(AmZ), m,m) + 


(m^ - A^)C'o(Ami^, a, m, m) - 2/3(3/2 Am - A)], 


Cii(A,to) = Cii(-Am'u, A,m,m) = ^ ^ [-Bo(-Am + A,m) + D{Am/ 2 + A)BoiAmV,m,m) 


{m^ - DA^)Co{Amv, a, m, m )], 


which in > 4 gives 


C'ii(A,m) = Cii{-AMV,A,m,m) = -[-Bo{-Am + A,m) + 4 (Am/2 + A)Bo{A^v,m,m)- 


{m^ - 4 :A‘^)Co{Amv, a, to, to) + 2/3(3/2Am - A)]. 


The calculation of the integral, 


/..A ^ 

Co(p, A, TOi,TO2) = -- 


d^k 


iTT^ J [kP- — nil + *^)((^ ~ pY ~ ^2 + • fc — A + ie) 

is done in [63]. For some calculations, we used the program FeynCalc [64]. 


Appendix B: Loop functions 


Loop functions entering Eq. (20) are listed here: 

1 


Ai{mi) = ((-5oo(Ap.,TOi) + Bii(Ap.,TOi) - AM.Bi(Ap.,TOi) 


—.Boo(Ap. + Am, mi) — i?ii(Ap). + Am, mi) — AmBi{Ad- + AM,mi))/2 


(5oo(Ap., Ap. + Am, mi) + B[^{Ap>, Ao* + Am, mi)) + 
+ gg {B'ooiAp, Ao + Am, mi) + 2i?Qg(Ap», Ap» + Am, mi))) , 


(mi) — i^ooi^P ”j + Am, mi) + Bh^Ap* , Ap» + Am, mi)) + 


+ gg {Booi^Pi^D + Am, mi) + 2Bqq(Ap., Ap. + Am, mi))) , 
for the Zlsi(2460)+ —>■ £)*+7r° decay mode and 

1 

IdiT^f 


A2{m,) = ^ ((Boo(Ap,TOj) + Sii(Ap,TOi) - AM.Bi(Ap,TOi) 


—.Boo(Ap + Am, mi) — Bii{Ap + Am, mi) — AmBi{Ad + Am, mi))l2 


(Bl) 


(B2) 
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the case of 

Ap [GeV] A 

p. [GeV] A 

D [GeV] Ap 

. [GeV] A 

M [Ge 

Al 

-0.59 

-0.45 

-0.06 

-0.04 

-0.35 

-4.2 

-0.48 

-0.31 

0.09 

O.l 

-0.33 

A'l 

-0.47 

-0.35 

-0.14 

0 

-0.35 


-0.33 

-0.21 

0 

-0.14 

-0.33 


Table I. Mass differences 


—h^ (^Bqq{Ap, Ap) + Am, nii) + B[i{Ap, Ap + Am, nii)) + 

+ ^ggBooi^D*,^p* + ^M,'mi)) , (B3) 


A2{m^) = 




{—h? (^Bqq{Ap, Ap + Am, mi) + B[i{Ap, Ap + Am, mi)) + 


+ iggBQQ{Ap-,Ap, + AM,mi)) , 


(B4) 


for the Z1 *q(2317)+ —-D+7 r° decay mode. 

The functions Bj which are present in Eq.(22) are: 

Bp+{m{) = {^g^BQoi'mp* - mjj+,Tni) - h'^{B'^Q{mp - m^+,TOi) + B[j^{mp - , 


Bjj*+[mi) = ^Q^2f2 {^9^i ~ Booimp - mjj*+,mi) + 2BQo(mD. - m^.+ ,TOi)) 


- h^{BQt^{mp> - mjj,+ ,mi) + B[^{mp> - 


B, 


= lQ^2f2 {^g^B'^Q{mp„ - mp^,mi) - h^[B'QQ{mp - mp^,m,) + B[^{mp - mp^ 


1 I I 5 


Bp: 


(mi) = ^Q^2f2 {^9^i ~ B'ooimp - mp:,mi) + 2BQo(mp. - mp.,mi)) + 


(B5) 


(B6) 

(B7) 


- /i^(i?oo(™D* - mp:,m{) + - mp:,mi))j . (B8) 

The expressions for Bj{mi) can be obtained from the above expressions of Bj{mi) by substituting masses of D mesons 
by the masses of Dg mesons. 

Here, Hi, Hqo and Hu, Hgo and HJi are loop integrals defined in Appendix A. Mass differences Ap, Ap, Ap. Ap. 
are defined as a mass differences between the appropriate state and the initial state while Am is the mass difference 
between final and initial state. Therefore, mass differences entering different amplitudes are not the same. The values 
are given in Table I. 
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